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The nonlinear probability distribution function in models 
with local primordial non-Gaussianity 
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ABSTRACT 

We use the spherical evolution approximation to investigate nonlinear evolution from 
the non-Gaussian initial conditions characteristic of the local f„i model. We provide 
an analytic formula for the nonlinearly evolved probability distribution function of the 
dark matter which shows that the underdense tail of the nonlinear PDF in the /„; 
model should differ significantly from that for Gaussian initial conditions. Measure- 
ments of the underdense tail in numerical simulations may be affected by discreteness 
effects, and we use a Poisson counting model to describe this effect. Once this has been 
accounted for, our model is in good quantitative agreement with the simulations. 
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1 INTRODUCTION 

The most common inflation model (the single scalar field, slow-roll inflation) predicts the primordial perturbations to be 
approximately described by the Gaussian statistics. Detections of non-gaussianity can discriminate between different inflation 

models, so the study of various cosmological probes of primordial non -gaussianity has attracted much recent attention. 

For example, the CMB has been used to constrain non-gaussianity IjKomatsu et al .I2OO8I : lHikage et al .l2008l : lYadav fc Wandel3 
20081 : iMcEwen et aLl boOS'l. In addition, various aspect s of large scale structure have also beeii exainined as probes of non 



Lo Verde et al. 



20081: iDalal et ahlboosl : IMatarrese fc Verde 



2008 



Carbone et al 



. ^ „ ^- .Slos^ et al 

200^ Desiacgues et al.l2008h: higher ord er effect s on the galaxy power sp ectrum |McDonald^^200S ; Taruya et al. II2OO8I); and 
the galaxy bispectrum jScoccimarro et al.i,2004 : ,Sefusatti fc KomatsiiboOT ). And. lsiosar. (,2008. ) has discussed how to optimize 
the constraint on /„; by applying weightings on subsamples from a single tracer. 



gaussianity. These include halo abundances and bias jLucchin et al.l 1988 : Robinson fc Bakei 2000l: Matarrese et al 



20081: Ufshordi fc ToUevlboOSl: 



200c 



Recently, iGrossi et al 



2008) used N-body simulations to study the effect of primordial non-gaussianity on the nonlinear 
probability distribution function (PDF) of the dark matter field. They found that departures from the Gaussian prediction 
are strongest in underdense regions and on small scales. The main goal of the present study is to provide some analytic 
understanding of their results. We do this by studying the nonlinear evolution of the dark matter PDF. When the initial 
conditions are Gauss ian, then the nonlinear evolution can be modelled by the spherical and ellipsoidal collapse models 
I Lam fc ShetiilboOSal lbh. In these models, the evolution of the PDF depends on two ingredients: a model for the dynamical 
evolution from an initial state to a final one, and the appropriate average over the initial states. In particular, the spherical and 
ellipsoidal collapse models approximate the nonlinear dynamics as 1^1 or 3— >1 mappings from the initial to the final state 
that do not depend on the Gaussianity of the initial conditions. Primordial non-Gaussianity enters only because it determines 
the initial set of states. As a result, essentially all of the machinery developed for the Gaussian case can be simply carried 
over to the present study. 

We recap the definitions of the local non-gaussian /„; model in Section [2. II The effect of smoothing the initial distribution 
is described in Section [?!2l Section[2]3]describes the nonlinear PDF associated with the spherical evolution model and compares 
the results with the Gaussian case. Section 12.41 compares our predictions with measurements in simulations. We summarise 
our results in Section |3l 
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2 THE PDF OF THE DENSITY IN THE LOCAL NON-GAUSSIAN MODEL 



In the fni mo dels of current in terest, the initial density fluctuatio n fleld is only mildly perturbed from the Gaussian. There- 
fore, following Lo Verde et al. (_2008) and Desiacques et al. ( 20081 ) . we use the Edgeworth expansion as a convenient way to 
summarize our results. Becau se we are interested in mass scales which are s ubsta ntially larger than that of a single collapsed 
halo - the regimes studied by Lo Verde et al. 1 20081 ) and Desiacques et al. ( 20081 ) - we are certainly in the regime where the 
Edgeworth expansion is a useful approximation. 



2.1 The unsmoothed initial distribution 

The primordial perturbation potential <1> of the local non-gaussian field is 

<E' = <^ + /n!(<^'-{<^')), (1) 

where </> is a Gaussian potential field and /„; is a scalar. We will use P(j>{k) to represent the power spectrum of <j); in what 
follows we will set P,j,{k) = AA:"""*, where ns ~ 1, and A is a normalization constant that is fixed by requiring that the rms 
fiuctuations in the associated non-Gaussian initial density field (which we will define shortly) have value as- 

We define M as the real, symmetric 3x3 tensor whose components are proportional to the second order derivatives of 
the potential <I>: 

t&ij = <l)ij + 2fnl{(l>i(t)j + 4>4>ij), (2) 

where = dicj) and — didjcj). We will sometimes refer to M as the shear or deformation tensor associated with the 
potential $. 

Correlations between the will be very useful in what follows. These depend on the correlations between (j> and its 
derivatives. However, because 4> is Gaussian, they can be computed easily. Specifically, 



^2 

{<t)xj4>) = ^^^^i {(pij(pk) = 0, {(l>ij(t)ki) = j^Sijki (4) 

where Sij is the dirac delta function, and we have defined 

Sijki = SijSki + SikSji + SilSjk, (5) 
cr| = J dkP^{k)]Vf{k)W^{kR)(-k'^y, (6) 

where 

3c2 



and T{k) is the CDM transfer function. (Strictly speaking, we are currently interested in the limit in which W — 1. We have 
defined the more general expression so that it can be used in the following subsections.) Thus, 

{$,,$^0 = ^5,,^, and (8) 

In 3 

{^ij ^kl $mn) = 2/„l-^— [SijSklmn + SklSijrnn + SmnSijkl] + 0(f„l), (9) 
Loo 



where 



„2 „2 

'^NG _ ^ I 4 j-2 

15 " 15 



„2 \ 2 „2„2 
15 



2.2 The smoothed initial field 



In the spherical evolution model, the quantity of interest is the initial overdensity Si smoothed on some scale V . Most of the 
complication in /„; models arises from the fact that the effect of smoothing is non-trivial. This nontriviality is a consequence 
of the fact that a smoothed Gaussian field is itself Gaussian, but this self-similarity does not hold for generic random fields. 

We now calculate the distribution function of the initial overdensity 5i smoothed on scale R in the fni model. Our goal 
is to approximate this initial PDF using the Edgeworth expansion: 



p{5i\R) ASi = 

where z^(-R) ~ Si / (jng{R) and H:i{v) — u{v^ 



2tt 



1 + 



GNG{R)S-i{R) 
6 



H,{u{R)) 



Au{R), 



(11) 



3). Therefore, we must specify how to calculate (tjvg and (tjvg 
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The Fourier transform of the initial overdensity is related to the Fourier transform of <1> by 
symmetry implies that the power spectrum and bispectrum of 5i are 



-k^ M{k) So spherical 



(27r)3 



Pi,(fc) = {-k^fM^(k)Pi,{k) = k^M'\k) 

Bs, (fcl, fc2, kl2) = {-kl)[-kl)[-kl2)M(k^)M(k2)M{kl2)B.i>{kl,k2, fcl2) 

B*(fci,fc2,fci2) = 2Ui [P0(fci)P^(fc2) + cyclic] 
From equations (|8} and ((9]), we have 

,2. 1 /■ dfc 



dqP4g)P^(|fc-q|) 



2 



2fnnNG{R) 



(Sf) 



(27r)3 y A; 

(5f) _ 2/„i7^G+0(/ni) 



(52)3/2 

2 

~(2^ 



47rfc'^M^ (A:) Rj (fc) (fci?) , 
where 



' NG 



I ^^i^C"^^) j ^klW{k2R) j A^ll2W{kl2R) Bs,{ki 



k2, ki 



(12) 

(13) 
(14) 



(15) 
(16) 
(17) 



with /ii2 = cos 012, and di2 is the angle between fci and fe2: ki2 = fci + ^2 + 2fcifc2 ^12- 

To second order in the variance <j%c{R) is the sum of two terms; one is the same as for the Gaussian, and the 
second, which is proportional to /^;, involves a convolution of the power spectrum with itself. This term has two infra-red 
singularities (at q = Q and k = q respectively) for tIs < 4. Fortunately, these can be removed by rewriting the power spectrum 
as 



Pi,{k) = P^k) + 



(27r)3 



dq[P4q)P^{\k-q 



P^{k)P^{q)-P^{k)P^{\k-q\)], (18) 
< 100, this renormalization is just a few percent effect, but 



where P4,{k) has a new normalization (|McDonaldll2008r ) . For 
the procedure is essential for removing the singularities. 

Equation (|17|l indicates that {S'f) scales li nearly with f^i- On sca les larger than about 100/i~^Mpc, the integral which 
defines 'yNG can be approximated analytically ( Scoccimarro et al.ll2004l ). but on smaller scales, the integral must be evaluated 
numerically. For models of interest, aSz is only weakly scale depende nt: e.g., it is ~ — 0.02(/n, i/100) on ~ 100^~^Mpc, and is 
less than a factor of two larger on scales (~ l/i~^Mpc) (e.g.. Figure 1 in lScoccimarro et al ] |2004h . This means that \cjng Sal <C 1 
on the scales of interest in this paper, justifying our use of the Edgeworth expansion. In addition, note that if aNa{R)S3{R) 
were independent of R, then the Edgeworth expansion would be a function of u only. In this case, the initial (non-Gaussian) 
PDF would be scale-independent, in the sense that the PDF would have the same functional form for all smoothing scales, 
just as it does for the Gaussian. Over a sufficiently narrow range of scales, this is a reasonable approximation. 



2.3 The PDF of the smoothed, evolved, nonlinear overdensity 



The previous section showed how to calculate the distribution of the initial overdensity 5i. Nonlinear evolution changes this 
distributions; this section shows how to estimate the evolved, nonlinear PDF of the overdensity. As noted in the introduction, 
this can be done by following the same steps as in Lam & Sheth (2008), but with the non-Gaussian initial distributions. For 
brevity, we show results for the spherical evolution model only. 

In the spherical model, the nonlinear overdensity in a region of volume V containing mass M is solely determined by the 
linear overdensity Si through a 1 1 mapping: 



pV 



'-I 



(19) 



( Bernardeau 1994 : Sheth 19981 ). where 5c ~ 5/3 and the exact value depends weakly on the cosmology. In this study we will 
use (5c = 1.66 which corresponds to the ACDM cosmology. The PDF of the nonlinear overdensity in a volume V associated 
with spherical collapse is 

din 1/ 



p'p{p\v) = PNG [5i (p) mp))^ ^^^p, 

where Png{5i\Vi) is the initial PDF of 5i at a given smoothing scale Vi. This with equation (|lip implies that 



(20) 



p%{,p\V) = 



P^Pg(p|V) 



exp 



H- 



Sfjp) 
2oHp) 

o{p)S3{p) 



Ha 



Slip) 



<^(P)S3(P) 



H3 



Slip) 

oip) 



oip) 



(21) 



where po denotes the smoothed nonlinear PDF assoociated with Gaussian initial conditions, (jip) — ctjvg(P = (3M/47rp)^'''') 
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Figure 1. Nonlinear overdensity PDF using the spherical collapse model in cells of radius 8/i^^Mpc. Long-dashed, solid, and dotted curves 
in the upper panel show ln(pp(p)) (from equation I21|l as a function of Inp for /„; = —100,0 and 100. Symbols show the PDF measured 
in the /„; = 100 simulation, and long-dashed-short-dashed curve shows the associated Poisson-sampled prediction (equation [22j . The 
lower panels show the log of the ratio between the /„; 7^ predictions and that for Gaussian initial conditions, for which /„; = 0. Filled 
and empty symbols are similar ratios of the measured PDFs for /„; = —100 and 100 respectively. 



and 7o- = — 3dlncr^/dlnM is not to be confused with the skewness parameter we defined earher. F inaUy, we set p' = Np and 
p'^p'ip') = P{p) to ensure that J dp'p'p'(p') and J dp' p'{p') both equal unity (see discussion in lLam fc Shethll2008b. .al . 

The non-gaussian modification contributes the final term of the right hand side of equation (|21|l . Since Si = when p = 1, 
there is little or no correction to the Gaussian case at p ~ 1. However, there is an effect at the low and high density tails. To 
see what it is, note that aSs is only a weak function of R, so it is a weak function of p, and hence the main p dependence is 
due to H3. To see what this dependence is, suppose that 70- — —6/5 (this is close to its actual value on the scales we consider 
when comparing with simulations in the next section). Then Si/a ~ (5/3)(l — p^'^^^)/ (ayp^^^^) ~ (5/3crv)(p^''^ — p~'^^^), 
where cry is the variance of the initial field when p = 1. When p ^ 1, then Si/a ~ (5/3(tv) p^^^, so H3 oc p'^^^. As a result, 
for /„; > 0, the high density tail of the PDF is increasing suppressed compared to the Gaussian case as p increases. At low 
densities, p ^ 1, then Si/a ~ — (5/3(tv) p~^^^, so H3 oc —p~^^^: the low density tail is enhanced for f„i > 0. The dependence 
on p is stronger for p ^ 1 than for p ^ 1, so we expect the underdense tail to be a good probe of /„(. For sufficiently large 
\Si/a\, the non-Gaussian piece can be negative, signaling that our truncation of the Edgeworth expansion was inappropriate. 
Fortunately, we are generally only interested in scales that are not significantly affected by this truncation problem. 
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Figure 2. Same as previous figure, but for cells of radius 4/i ^Mpc. 



2.4 Comparison with simulations 

We compare the predi ctions of our approach with measurements of the nonhnear PDF in numerical simulations from 
Desiacaues et al. ( 20081 ) . These followed the evolution of 1024"^ particles in a periodic cube of sides 1600 /i~^Mpc. The back- 
ground cosmology is spatially flat, dominated by a cosmological constant, having {flrn,iib,na,h,as) = (0.279, 0.0462, ris = 
0.96,0.7,0.81), so the particle mass was 3 x 10^^h~^ Mq. The simulations sample the density field using discrete particles. 
This produces discreteness effects which are largely irrelevant, except in the least dense tails of the PDF. (E.g., the average 
number of particles in spherical cells with radius 4/i~^Mpc is ~ 60. So discreteness effects are severe at p < 1/6.) We account 
for this using the Poisson model: 

(22) 



p(iV|V) = / dM p{M\V) p{N\M) 



where p(N\M,V) = {M/rUp)^ exp{-M/mp)/N\, with equal to the particle mass (|Shethlll996l '). We then set p = N/nV 



when plotting the results. 

The symbols in the top panel of Figure[l]show the measured PDF for counts in spheres of radius 8/!.~^Mpc for /„; — 100. 
The dashed curve which is closest to these symbols shows equation (I22p for /„; = 100. The dashed, solid and dotted curves 
show the predictions of equation pi|l for /„; = —100, and 100. The bottom panel shows the ratio of the counts in the /„( 
models to those when /„; = 0. The outer set of curves show equation pif) and the inner set show the effect of accounting 
for discreteness effects using equation (|22p . This shows that the PDF for positive /„; is slightly skewed towards underdense 
regions, the reverse is true for negative /„;, and our model provides an excellent description of these trends. 

Figure [2] shows a similar analysis of the PDF on smaller scales, for which the discreteness effects are more pronounced: 
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at small p, the symbols in the top panel lie well above the curves associated with equation (|2H) . However, our Poisson model 
for the discreteness effect appears to be rather good: the prediction associated with equation ((22} provides a good description 
of the measurements. The bottom panel shows that discreteness effects tend to wash out the differences between the different 
f„i runs, but that our Poisson model does an excellent job of accounting for this effect (the predictions at small p are almost 
indistinguishable from the measurements). 



3 DISCUSSION 

We used the spherical evolution model to study the nonlinear evolved probability distribution function of the dark matter 
density field in the local primordial non-Gaussian /„; model. (The spherical model is able to provide a good description 
of the nonlinear PDF when /„; — 0.) For currently acceptable values of f„i, our approach shows that the signatures of 
primordial non-Gaussianity are small, but are most evident in underdense regions (equation [2T] and related discussion). The 
PDF measured in simulations can be affected by discreteness effects, especially in small underdense cells. The effect of this 
can be approximated by using a Poisson counting model (eauation l22p . Once this is done, our model is in very good agreement 
with measurements in numerical simulations (Figures [T] and [J]! , so we hope our equation (|2ip will be useful in studies which 
require knowledge of the nonlinear evolved PDF. 

Recently, following the same logic for why cluster abundances should be good probes of primordian non-Gaussianity, 
Kamionkowski et al. ( 20081 ) have suggested that void abundances should also be good probes. Our results provide further 



motivation for st udying underdense regions. We are in the process of developing a more complete model of voids and void 



shapes f following ISheth fc van de Wevgaertll2004l ') 



Finally, note that the parameter which describes the non-Gaussianity in the smoothed initial and final fields, a S3, is only 
weakly scale-dependent. Therefore, our analysis indicates that the non-gaussian distribution of the resulting nonlinear PDF 
(equation [21} can be written in terms of the scaled variable f = Si /a. So one could formulate a reconstruction of the initial 



fni field analogously to how this is done for the Gaussian case l|Lam fc Shethll2008al ). We have not pursued this further. 



ACKNOWLEDGEMENTS 

We thank V. Desjacques for help with measuring the PDFs in his simulations and R. Scoccimarro and E. Komatsu for many 
helpful discussions. Thanks also to S. Cole, K. Dolag, M. Grossi, W. Hu, Y. P. Jing, Z. Ma, T. Nishimichi, and R. Scoccimarro 
for discussions about resolution effects on the setting up of initial conditions in simulations. 



REFERENCES 

Afshordi N., ToUey A. J., 2008, ArXiv e-prints, astro-ph/0806.1046 

Bernardeau F., 1994, Astron. Astrophys., 291, 697 

Cole S., Kaiser N., 1989, Mon. Not. R. Astron. Soc, 237, 1127 

Carbone C, Verde L., Matarrese S., 2008, Astrophys. J. Lett., 684, LI 

Dalai N., Dore O., Huterer D., Shirokov A., 2008, Phys. Rev. D, 77, 123514 

Desjacques V., Seljak U., Iliev 1., 2008, Mon. Not. R. Astron. Soc, submitted (arXiv:0811.2748) 
Grinstein B., Wise M., 1986, Astrophys. J., 310, 19 

Grossi M., Branchini E., Dolag K., Matarrese S., Moscardini L., 2008, ArXiv e-prints, astro-ph/0805.0276 

Hikage C, Matsubara T., Coles P., Liguori M., Hansen F. K., Matarrese S., 2008, Mon. Not. R. Astron. Soc, 389, 1439 

Kaiser N., 1984, Astrophys. J. Lett., 284, L9 

Kamionkowski M., Verde L., Jimenez R., 2008, ArXiv e-prints, astro-ph/0809.0506 

Komatsu E., Dunkley J., Nolta M. R., Bennett C. L., Gold B., Hinshaw G., Jarosik N., Larson D., Limon M., Page L., 
Spergel D. N., Halpern M., Hill R. S., Kogut A., Meyer S. S., Tucker G. S., Weiland J. L., WoUack E., Wright E. L., 2008, 
ArXiv e-prints, astro-ph/0803.0547 

Lam T. Y., Sheth R. K., 2008a, Mon. Not. R. Astron. Soc, 386, 407 

Lam T. Y., Sheth R. K., 2008b, Mon. Not. R. Astron. Soc, 389, 1249 

Lo Verde M., Miller A., Shandera S., Verde L., 2008, Journal of Cosmology and Astro-Particle Physics, 4, 14 
Lucchin F., Matarrese S., Vittorio N., 1988, Astrophys. J. Lett., 330, L21 
Matarrese S., Verde L., 2008, Astrophys. J. Lett., 677, L77 
Matarrese S., Verde L., Jimenez R., 2000, Astrophys. J., 541, 10 
McDonald P., 2008, ArXiv e-prints, astro-ph/0806.1061 

McEwen J. D., Hobson M. P., Lasenby A. N., Mortlock D. J., 2008, Mon. Not. R. Astron. Soc, 388, 659 



Dark matter PDF in the fni model 



Mo H. J., White S. D. M., 1996, Mon. Not. R. Astron. Soc, 282, 347 

Robinson J., Baker J., 2000, Mon. Not. R. Astron. Soc, 311, 781 

Scoccimarro R., Sofusatti E., Zaldarriaga M., 2004, Phys. Rev. D, 69, 103513 

Sefusatti E., Komatsu E., 2007, Phys. Rev. D, 76, 083004 

Sheth R. K., 1996, Mon. Not. R. Astron. Soc, 279, 1310 

Shcth R. K., 1998, Mon. Not. R. Astron. Soc, 300, 1057 

Sheth R. K., Tormen G., 1999, Mon. Not. R. Astron. Soc, 308, 119 

Sheth R. K., van de Weygaert R., 2004, Mon. Not. R. Astron. Soc, 350, 517 

Slosar A., 2008, ArXiv c-prints, astro-ph/0808.0044 

Slosar A., Hirata C, Seljak U., Ho S., Padmanabhan N., 2008, Journal of Cosmology and Astro-Particle Physics, 8, 31 
Taruya A., Koyama K., Matsubara T., 2008, ArXiv e-prints, astro-ph/0808.4085 
Yadav A. P. S., Wandelt B. D., 2008, Physical Review Letters, 100, 181301 



